Neutron Star Structure and the Neutron Radius of 208 Pb 
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We study relationships between the neutron-rich skin of a heavy nucleus and the properties of 
neutron-star crusts. Relativistic effective field theories with a thicker neutron skin in 208 Pb have a 
larger electron fraction and a lower liquid-to-solid transition density for neutron-rich matter. These 
properties are determined by the density dependence of the symmetry energy which we vary by 
adding nonlinear couplings between isoscalar and isovector mesons. An accurate measurement of the 
neutron radius in 208 Pb — via parity violating electron scattering — may have important implications 
for the structure of the crust of neutron stars. 



It is an extrapolation of 18 orders of magnitude from 
the neutron radius of a heavy nucleus — such as 208 Pb 
with R n w 5.5 fm — to the approximately 10 km radius 
of a neutron star. Yet both radii depend on our incom- 
plete knowledge of the equation of state of neutron-rich 
matter. Therefore, an accurate measurement of the neu- 
tron radius in 208 Pb may have important implications for 
some neutron star properties. 

Heavy nuclei are expected to have a neutron-rich 
skin. This important feature of nuclear structure arises 
because of the large neutron excess and because the 
Coulomb barrier reduces the proton density at the sur- 
face. The thickness of the neutron skin depends on the 
pressure of neutron-rich matter: the greater the pressure, 
the thicker the skin as neutrons are pushed out against 
surface tension. The same pressure supports a neutron 
star against gravity Q. Thus models with thicker neu- 
tron skins often produce neutron stars with larger radii. 

Neutron stars are expected to have a solid crust of 
nonuniform neutron-rich matter above a liquid mantle. 
The phase transition from solid to liquid depends on the 
properties of neutron-rich matter. Indeed, a high pres- 
sure implies a rapid rise of the energy with density mak- 
ing it energetically unfavorable to separate uniform mat- 
ter into regions of high and low densities. Thus a high 
pressure typically implies a low transition density from a 
solid crust to a liquid mantle. This suggests an inverse 
relationship: the thicker the neutron-rich skin of a heavy 
nucleus, the thinner the solid crust of a neutron star. 

In this letter we study possible "data-to-data" relations 
between the neutron-rich skin of a heavy nucleus and the 
crust of a neutron star. These relations may impact neu- 
tron star observables. Indeed, properties of the crust are 
important for models of glitches in the rotational period 
of pulsars (||| > for the shape and gravitational radiation 
of non-spherical rotating stars Q and for neutron-star 
cooling ||. Note that the skin of a heavy nucleus and 
the crust of a neutron star are composed of the same 
material: neutron-rich matter at similar densities. 



The Parity Radius Experiment (PREX) at the Jeffer- 
son Laboratory aims to measure the neutron radius in 
208 Pb via parity violating electron scattering Hf7|. Par- 
ity violation is sensitive to the neutron density because 
the Z° boson couples primarily to neutrons. The result 
of this purely electroweak experiment could be both ac- 
curate and model independent. In contrast, all previous 
measurements of bulk neutron densities used hadronic 
probes that suffer from controversial uncertainties in the 
reaction mechanism (see for example Ref. ||). PREX 
should provide a unique observational constraint on the 
thickness of the neutron skin in a heavy nucleus. In this 
letter we explore some of the implications of this mea- 
surement on the structure of neutron stars. 

Microscopic calculations of the energy of neutron mat- 
ter constrain both the neutron skin in 208 Pb and the crust 
of a neutron star; see for example Ref. ||. However, 
these calculations of infinite neutron matter are not di- 
rectly tested by observable properties of finite nuclei such 
as their charge densities or binding energies. Moreover, 
nonrelativistic calculations of symmetric nuclear matter 
have not succeeded in predicting the saturation density. 
It thus becomes necessary to fit some properties of a 
three-body force in order to reproduce nuclear saturation. 
Indeed, the properties of A = 8 pure neutron drops calcu- 
lated in Ref. |Hj may depend on the three-nucleon force 
used. Thus we feel that it is important to distinguish di- 
rect finite-nucleus measurements — such as PREX — from 
theoretical neutron-matter "observables" based solely on 
calculations. Indeed, PREX may provide an important 
test of these calculations Jll]] . 

We start with a relativistic effective field theory |l^ | 
that provides a simple description of finite nuclei and 
a Lorentz covariant extrapolation for the equation of 
state of dense neutron-rich matter. The theory has an 
isoscalar-scalar 4> (sigma) meson field and three vector 
fields: an isoscalar V (omega), an isovector b (rho), 
and the photon A. We now supplement the Lagrangian 
with new nonlinear sigma-rho and omega-rho couplings. 
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These couplings allow us to change the density depen- 
dence of the symmetry energy which changes both the 
thickness of the neutron skin in 208 Pb and the neutron- 
star crust. 

The interacting Lagrangian density is given by p2[ , 
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We consider a number of different parameter sets. 
First, we note that the nonlinear rho coupling £ will mod- 
ify the density dependence of the rho mean field and this 
could change the neutron-skin thickness. However, un- 
less £ is made very large, this term was found to have 
a small effect |T^| . Thus for simplicity, we set £ = 
in all our parameter sets. Second, note that we could 
have added a cubic sigma-rho interaction of the form: 
£ 3 = MK s {g s (f)){g' 2 p h ll • b p ). A nonzero A 3 does change 
the thickness of the neutron skin in 208 Pb — but at the 
expense of a change in the proton density. Therefore, we 
set A 3 = and focus exclusively on A4 and A v . 

We start with the original NL3 parameter set of 
Lalazissis, Konig, and Ring Jl3[] . (Note that a small ad- 
justment of the NNp coupling constant was needed to 
fit the symmetry energy of nuclear matter at a Fermi 
momentum of fe^ = 1.15fm~ ; see text below). The NL3 
set has C = A4 = A v = and provides a good fit to the 
ground-state properties of many nuclei. In this model 
symmetric nuclear matter saturates at kp — 1.30 fm 
with a binding energy per nucleon of £?A4= —16.25 MeV 
and an incompressibility of K = 271 MeV. All other pa- 
rameter sets considered here have been fixed to the same 
saturation properties. 

We now add the new nonlinear couplings A4 and/or A v 
between the isoscalar and the isovector mesons. At the 
same time we adjust the strength of the NNp coupling 
constant (g p ) from its NL3 value to maintain the symme- 
try energy of nuclear matter unchanged (see text below). 
Note that neither A4 nor A v affect the properties of sym- 
metric nuclear matter since b M = 0. Hence, the saturation 
properties remain unchanged. Our goal is to change the 
neutron density and the neutron-skin thickness in 208 Pb 
while making very small changes to the proton density 
which is well constrained by the measured charge den- 
sity |Q. The two new couplings (A 4 and A v ) change 
the skin thickness in 208 Pb by similar amounts. Yet they 
have different high-density limits [Q. For A v = the 
symmetry energy is proportional to the baryon density 
p in the limit of very high density, while it only grows 
like p 1 / 3 for nonzero A v . This can produce noticeable 
differences in neutron star radii, as we show below. 

Let us start with A4 = 0. For a given omega- rho cou- 
pling A v we readjust only the NNp coupling constant 



g p in order to keep an average symmetry energy fixed. 
The symmetry energy at saturation density is not well 
constrained by the binding energy of nuclei. However, 
some average of the symmetry energy at full density and 
the surface energy is constrained by binding energies. As 
a simple approximation we evaluate the symmetry en- 
ergy not at full density kp ~ 1.30 fm -1 but rather at 
an average density corresponding to kp = 1.15 fm -1 . 
Thus, all our parameter sets have a symmetry energy 
of 23.50 MeV at kp = 1.15 fm -1 . Note that the original 
NL3 parameter set predicts a symmetry energy of 37.4 
MeV at full saturation density and close to 23.50 MeV at 
kp = 1.15 feu -1 fl3|| . This simple procedure produces a 
nearly constant binding energy per nucleon for 208 Pb as 
A v is changed, as can be seen in Table |. Moreover, Ta- 
ble I shows that increasing A v reduces the neutron-skin 
thickness significantly — while maintaining the proton ra- 
dius nearly constant. In the following we plot our results 
for a range of A v values for which the proton radius is 
within 0.01 fm of its A v = value. 

To study the solid crust of a neutron star we make a 
simple random-phase-approximation (RPA) calculation 
of the transition density below which uniform neutron- 
rich matter becomes unstable against small amplitude 
density fluctuations. This provides a lower bound to the 
true transition density [Gl|. We start with the longi- 
tudinal dielectric function e^, as defined in Eq. (68) of 
Rcf . |ll| , evaluated at an energy transfer qo = and at 
an arbitrary momentum transfer q. That is, 



e L {qo = 0,q) = det (1 - D L U L ) 



(2) 



Here is a longitudinal polarization matrix describing 
particle- hole excitations of a uniform system of protons, 
neutrons, and electrons in beta equilibrium, as given in 
Eq. (56) of Ref. |lq ]. The matrix Dl, describing me- 
son and photon propagation, follows from Eq. (57) of 
Ref. jl6) — but includes additional terms to account for 
the nonlinear nature of the meson self-couplings [fj7| . 
We estimate the transition density p c by computing the 
largest density at which 6^(0, q) <0 for any given q. 

TABLE I. Results for the NL3 Parameter Set. The bind- 
ing energy per particle in 208 Pb is B.E., R p is the proton and 
Rn — R P is the difference between neutron and proton radii 
in Pb. (Note that we do not include center of mass correc- 
tions.) Finally, p c is our estimate for the transition density of 
neutron-rich matter from a nonuniform to uniform phase. 
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FIG. 1. Estimate of the transition density from nonuni- 
form to uniform neutron-rich matter versus neutron- 
minus-proton radius in 208 Pb. The curves are for the 
four parameter sets described in the text. 

In Fig. 1 we display the transition density for various 
parameter sets (see Table ||) as a function of the pre- 
dicted difference in the root-mean-square neutron and 
proton radii Rn — Rp in 208 Pb. The curves are param- 
eterized by different values of A v , as shown in Table Q. 
The NL3 parameter set saturates nuclear matter with 
a relatively small value of the nucleon effective mass: 
M* = M-g s <j> = 0.59M. The parameter set S271 satu- 
rates nuclear matter as NL3 but with M* = 0.70M. This 
set also has C = 0. The two remaining curves in the figure 
are for parameter sets having £ = 0.06 and both satu- 
rate nuclear matter with M * = 0.80M. (Set Z271v has 
a nonzero A v while set Z2714 uses a nonzero A4). Note 
that the scalar mass m s for parameter sets S271, Z271v, 
and Z2714 is adjusted to reproduce the proton radius in 
208 Pb as computed with NL3. Figure 1 displays a clear 
inverse correlation between the transition density and the 
neutron-skin thickness R n — R p . The transition density 
expressed in fm~ 3 is about, 



Pc 



0.16 - 0.39 (Rn - Rp) , 



(3) 



with the skin thickness expressed in fm. Moreover, this 
correlation seems to be insensitive to M* or to using A4 
or A v to change R n — R p . These results suggest that a 
measurement of the neutron radius in 208 Pb will provide 
considerable information on the transition density. 

TABLE II. Model parameters used in the calculations. 
The parameter k and the scalar (m a ) and vector (m v ) masses 
are given in MeV. The nucleon and rho masses are kept fixed 
at M — 939 and m p = 763 MeV, respectively. 
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6.68 
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6.17 


0.15634 
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Note that Fig. 1 only shows our results. Yet all 
other calculations that we are aware of also give con- 
sistent results. For example, the nonrelativistic micro- 
scopic equation of state of Friedman and Pandharipande 
has a transition density of p c = 0.096 fm~ 3 according to 
Lorenz et al. M. For this equation of state Brown finds 
Rn - -Rp = 0.16 ± 0.02 fm Q. These numbers are in 
excellent agreement with Eq. (3). 

Brown also finds a linear relation between R n — R p 
and the derivative of the energy of neutron matter versus 
density dE/dp evaluated at p = 0.1 fm -3 [^T). He consid- 
ers a large variety of nonrelativistic Skyrme interactions. 
Our results for dE/dp versus R n — R p are completely 
consistent. We expect these common dE/dp values to 
give similar p c values consistent with Eq. (3) for these 
Skyrme interactions. This is because dE/dp is related to 
the pressure while p c depends on the density dependence 
of the pressure. 

Finally, for the relativistic interaction TM1 of Suga- 
hara and Toki fl~8| , we calculate from Eq. (2) p c m 0.059 
fm" 3 and R n - R p = 0.27 fm. Again, the numbers are 
again in good agreement with Eq. (3). 

In Fig. 2 we show the electron fraction per baryon Y e 
versus density for uniform neutron-rich matter in beta 
equilibrium. We include results only for the S271 pa- 
rameter set as all other sets yield similar results. The 
different curves are for different values of A v which pre- 
dict the indicated R n — R p values. The curves start near 
the transition densities displayed in Fig. 1. The elec- 
tron fraction Y e is determined by the symmetry energy 
while R n —R p is sensitive to the density dependence of the 
symmetry energy. Therefore a measurement of R n — R p 
constrains the growth of Y e with density. If R n — R p is 
greater than about 0.24 fm, Y e becomes large enough to 
allow the direct URCA process Jl9| to cool down a 1.4 
solar mass neutron star. 
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FIG. 2. Electron fraction Y e versus baryon density for 
uniform neutron-rich matter in beta equilibrium using 
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the S271 parameter set (other sets yield similar results). 
The curves are for different values of A v that predict the 
indicated values of R n — Ftp for 208 Pb, in fm. 

We now consider the radius R of a 1.4 solar mass neu- 
tron star. For a given parameter set, R increases with 
R n — R p . However, as one changes the parameter set to 
increase M* or £ the equation of state becomes softer at 
high density. As a result R decreases for fixed R n — R p . 
Also, the high density equation of state is softer with A v 
than with A4 so Z271v gives slightly smaller stars than 
parameter set Z2714. We conclude that R is not uniquely 
constrained by a measurement of the neutron-skin thick- 
ness because R n — R p only depends on the equation of 
state at normal and lower densities while R is also sen- 
sitive to the equation of state at higher densities. Yet 
one may be able to combine separate measurements of 
R n — R p and R to obtain considerable information about 
the equation of state at low and high densities. For ex- 
ample, if R n — R p is relatively large while R is small this 
could suggest a phase transition. A large R n —R p implies 
that the low density equation of state is stiff while a small 
R suggests that the high density equation of state is soft. 
The transition from stiff to soft could be accompanied by 
a phase transition. 




a measurement of the skin thickness will constrain the 
density dependence of the symmetry energy. 

2) The density dependence of the symmetry energy is 
adjustable in our relativistic effective field models while 
still reproducing nuclear-matter properties and other 
ground-state observables. Indeed, our models can pro- 
vide a Lorentz-covariant extrapolation for the high den- 
sity equation of state with a symmetry energy that rises 
slower with density relative to earlier relativistic mean- 
field models. 

3) The electron fraction Y e of neutron-rich matter 
in beta equilibrium is correlated with the neutron-skin 
thickness in 208 Pb. The thicker the neutron skin the 
faster Y e rises with density. In our models a neutron- 
skin thickness of the order of 0.24 fm or larger suggests 
that Y e will become large enough to allow a direct URCA 
process to cool down a 1.4M© neutron star. 

4) We have found an inverse correlation between the 
neutron-skin thickness and the density of a phase tran- 
sition from nonuniform to uniform neutron-rich matter. 
In our models the transition density obeys the approxi- 
mate relation: p c « 0.16 — 0.39 (R n — R p ), with R n — R p 
in fm and p c in fm -3 . This suggests that a neutron skin 
measurement in 208 Pb can provide important informa- 
tion on the thickness and other properties of the crust of 
a neutron star. 

5) Microscopic calculations of the energy of neutron 
matter constrain the density dependence of the sym- 
metry energy and hence the neutron-skin thickness in 
208 Pb. Therefore a neutron skin measurement may pro- 
vide an important observational check on such calcula- 
tions that is not provided by other observables. More- 
over, a neutron-skin measurement may constrain three- 
body forces in neutron-rich matter. 

We acknowledge useful discussions with A. Brown, R. 
Furnstahl, J. Lattimer, S. Reddy and B. Serot. We thank 
A. Brown for sharing results prior to publication. This 
work was supported in part by DOE grants DE-FG02- 
87ER40365 and DE-FG05-92ER40750. 
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FIG. 3. Radius of a 1.4 solar mass neutron star ver- 
sus neutron-minus-proton radius in 208 Pb for the four 
parameter sets described in the text. 

In conclusion: 1) It is possible to fit nuclear 
observables — such as charge densities, binding energies, 
and single particle spectra — with effective field theories 
that predict a range of neutron-skin thicknesses. This can 
be done by adding nonlinear couplings between isoscalar 
and isovector meson fields or, in general, by adding inter- 
actions that modify the density dependence of the sym- 
metry energy. We conclude that the neutron-skin thick- 
ness is not tightly constrained by these observables. Yet 
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